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abstract 


The  reflection  and  scattering  of  acoustic  signals  by  the  bottom  of 
the  ocean  depend  upon  the  physical  properties  of  the  bottom  and  the  wave¬ 
length  spectrum  of  the  signals.  On  the  basis  of  the  frequency  and  gracing 
angle  dependence  of  the  reflected  or  scattered  signals,  physical  properties 
such  as  the  densities,  acoustic  velocities  of  subbottom  layers,  and  the 
rms  roughness' of  the  layers  can  be  estimated. 
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PREFACE 


This  report  has  been  written  in  two  parts.  Each  of  the  parts  is 
quite  independent  and  yet  each  part  covers  much  of  the  same  material. 
Their  difference  lies  in  that  Part  I  is  a  development  of  the  subject, 
whereas  in  Part  II  certain  simplified  results  are  applied  directly  to 
marine  geophysical  problems. 

While  Part  II  was  written  first  and  was  intended  to  assist  in  the 
interpretation  of  the  marine  geophysical  surveys,  the  existence  of  Part  I 
as  a  reference  was  assumed.  Part  I  formed  the  basis  of  an  invited  paper 
presented  at  the  74th  Meeting  of  the  Acoustical  Society  of  America,  Miami 
Beach,  Florida,  November  1967. 


C.  S. 


Clay 


t 


C.  S.  Clay,  "Measurement  of  the  ocean  bottom,  "  J.  Acoust.  Soc.  Am. 
42,  1183  (A)  (1967).  (HL  Inf.  Doc.  No.  139.) 
i  t 

Dr.  Clay  is  now  Professor  of  Geophysics,  Department  of  Geology 
and  Geophysics,  University  of  Wisconsin,  Madison,  Wisconsin. 
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PART  I.  THEORY  AND  EXPERIMENT 


1.  Interaction  of  Sound  at  the  Bottom  -  Interpretation 

Our  purpose  is  to  examine  rather  carefully  the  interpretation  of 
measurements  of  the  interaction  of  acoustic  signals  with  the  bottom  of  the 
ocean.  We  have  avoided  the  words  ''reflected,  "  "scattered,  "  or  "diffracted" 
because  these  words  describe  a  particular  mechanism  and  method  of  data 
analysis. 

The  basic  measurement  is  made  as  shown  in  Fig.  1.  The  source 
transmits  a  transient  signal  so  that  the  sound  that  interacts  with  the  bottom 
can  be  separated  from  other  paths  such  as  the  direct  path  or  multiple  sur¬ 
face  and  bottom  paths.  The  bottom  has  been  drawn  as  being  a  black  box, 
and  we  don1 1  know  much  about  it  except  that  it  is  passive,  It  is  also 
bounded  in  the  sense  that  we  can  confine  the  sound  field  to  a  given  region 
of  the  bottom. 


Fig.  1.  Interaction  of  sound 
with  the  ocean  bottom.  S  is 
the  source,  C  is  the  calibra¬ 
tion  receiver,  and  rev  is  the 
receiver.  The  ocean  bottom 
is  inside  an  imaginary  black 
box. 
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We  now  have  a  measurement  of  the  acoustical  signal  at  the  receiver. 


The  positions  of  the  source,  receiver,  and  insonified  region  are  known. 
Signal  processing  equipment  for  filtering,  correlating,  squaring,  displaying, 
etc.  is  available  and  so  we  need  not  discuss  that.  The  question  is  how  can 
we  use  the  geometry  of  the  measurement  to  measure  the  interaction  at  the 
bottom.  We  stress  the  word  "how"  because  the  way  we  chose  to  work  up 
the  data  requires  an  assumption  about  the  nature  of  the  interaction. 

The  procedure  we  follow  is  trial  and  error.  That  is,  we  assume 
a  description  of  the  bottom  and  then  calculate  theoretical  received  signals 
(with  or  without  signal  processing)  for  a  number  of  different  geometries. 
Depending  upon  how  good  the  comparison  of  the  theoretical  and  experimental 
measurements  is,  we  use  the  description  or  try  again. 

2.  Reflecting  Bottom 

Let  us  assume  the  box  to  be  a  reflector  for  the  first  trial.  The 

consequence  of  this  assumption  is  shown  in  Fig.  2.  By  means  of  the 

geometry  and  calibration,  we  compute  the  ratio  of  the  measured  signal 

p  to  the  theoretical  signal  p  ,  for  a  unit  perfect  reflector.  The 

meets  till  r 

ratio  is  termed  the  apparent  reflection  coefficient  ffit  .*  We  use  the 

ct 

word  "appareht"  because  B.a  is  the  result  of  an  interpretation  or  assump¬ 
tion,  and  not  a  direct  measurement. 


In  text  and  equations  ffi,  substitutes  for  script  R  . 
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Fig.  2.  Reflection  at  the  bottom 
pQ  is  the  signal  at  unit  distance 


Fig.  3.  Scattering  at  the  bottom; 
A  is  the  illuminated  area  and  its 
dimensions  are  much  less  than 
Rj  and  R^  . 


3.  Scattering  Bottom 

We  could  have  assumed  the  bottom  to  be  so  rough  that  it  scatters 
energy  into  all  directions.  The  scattering  interpretation  is  shown  in  Fig.  3. 
For  illustration,  the  density-velocity  contrast  at  the  bottom  is  assumed  to 
be  so  large  that  each  small  facet  of  the  bottom  can  be  regarded  as  being  a 
perfect  reflector.  The  apparent  scattering  coefficient  is  S  .  As  in 
the  previous  example,  S  is  the  result  of  an  interpretation.  Comparison 

ft 

of  the  dependence  of  p  .  and  p  .  upon  the  geometry  shows  that 

will  r  wii|  8 

JB-a  is  a  function  of  ^R^R^”*  whereas  is  a  function  of 

KV)  • 

These  two  very  different  interpretations  can  be  applied  to  the  same 
measurement.  It  is  obvious  that  we  need  to  know  whether  the  bottom  is 
smooth  or  rough  so  as  to  know  how  to  work  up  our  data.  Let  us  look  at 
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a  geological  and  geophysical  conception  of  the  bottom  of  the  ocean  -  as 
shown  in  Fig.  4.  This  somewhat  artistic  drawing  is  based  upon  a  great 
deal  of  knowledge.  Most  of  it  was  obtained  by  means  of  echo  sounders, 
bottom  samples,  deep  sea  cameras,  and  seismic  profilers.  Generally, 
it  is  a  good  picture,  but  we  probably  should  not  use  it  to  describe  the 
bottom  precisely  enough  to  predict  the  result  of  a  particular  measurement. 
The  reason  is  that  our  measurements  are  sensitive  to  details  that  are  .iot 
ordinarily  shown  or  measured.  These  charts  can  be  used  as  a  guide  for 
our  first  guess,  but  in  the  end  the  measurements  must  be  used  to  make 
a  description  of  the  bottom. 


TO*  65* 


Fig.  4.  Heezen  and  Tharp 
pictographic  chart  of  the  ocean 
bottom. 


4.  Reflection  and  Scattering  of  Sound  by  ah 
Irregular  Interface  -  Theory 

The  theory  of  scattering  of  waves  by  an  irregular  surface  has  been 

the  subject  of  a  number  of  papers  and  is  treated  in  several  books.  The 

2  3 

first  formulation  was  given  by  Rayleigh.  Brekhovskikh  and  Eckart 
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have  given  solutions  based  on  the  Helmholtz  equation.  Eckart' s  formula¬ 
tion  is  particularly  convenient  for  the  study  of  reflection  and  scattering  by 
a  randomly  rough  surface. 

4 

The  development  given  by  Tolstoy  and  Clay  is  essentially  that  of 
Eckart,  and  it  is  convenient  to  follow  their  notation.  We  make  the 
simplifying  assumption  that  the  source  and  receiver  are  very  far  away 
from  the  illuminated  area.  On  neglecting  the  slopes  of  the  surface,  their 
Eq.  (6.  19)  is 


ik(VRz) 

e  **  C  r»  n>  2i(ax+(5y+yr)  ,  , 

mar, —  J  DM'<e  pywdydx 


2a  = 

k  [sin  0^-sin  0^  cos 

2p  = 

-k  |sin  0^  sin  B^j 

2y  - 

-k  [cos  0  +cos  By  j 

K2  = 

•  V 

k  = 

Zv/\  =  w  /c 

B2  = 

SI  p  c(2ir)  1 

D  = 

illumination  function 

o  - 

■  m  ““ 

source  power 

Rj  is  in  the  x-z  plane 

6 ^  is  the  angle  between  the  x-z  plane  and 
plane  containing  and  the  z  axis 

£  is  the  elevation  relative  to  z=0 

JR.  is  the  reflection  coefficient 


(1) 


-6- 


Fig.  5.  Geometry  for 
Helmholtz  formulation. 


The  geometrical  quantities  are  defined  in  Fig.  5.  We  notice  that 
2a  ,  Zp  ,  Zy  are  the  differences  of  the  wave  numbers  along  Rj  and 
.  In  our  approximation,  they  are  constant  for  particular  source  and 
receiver  positions.  In  the  specular  direction,  the  horizontal  components 
(a,  |3)  are  zero  and  this  considerably  simplifies  the  integral: 


Themean  square  acoustic  pressure  is 


<PP*>  *  I  I  I  I  Da’  (e2i'll’=~C))  dy  dx  dy"  dx"  .  (3) 

4*  RjR2  '  ' 

The  average  is  taken  over  the  random  surface.  We  assume  the  distributions 
of  £  and  to  be  Gaussian  and  that  difference  (£' -4")  depends  upon 

their  separation.  Thus,  the  substitution  of  the  bivariate  Gaussian  function 
and  the  evaluation  of  the  average  yield  the  following: 
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\ 


with 

W  = 

1 

;2+c«2_2CC//iJi 

1  (4) 

2irc-2(l-+2)  6XP 

(ifcTj 

4/  = 

(5) 

<r  = 

rms  roughness. 

<pp*) 


y2jr2b2 

4’'2r1«2 


////  tireV'!M) 


dy  dx  dy"  dx" 


(6) 


We  have  chosen  to  use  a  Gaussian  type  of  illumination  function 
because  it  resembles  the  main  lobes  of  array  responses  and  has  no  side 
lobes.  With  the  aid  of  the  assumption  that  x  depends  only  upon  (x-x") 
and  (y-y  )  ,  the  following  substitutions  permit  the  integration  of  (6) 


D  =  e-(x2/2X2)-(y2/2Y2) 


x"  =  x-| 


y"  =  y-Ti 


D  -  e--2/*2 


(6a) 


<PP*> 

.  y  a2B2 

4itR2R 

or 

<PP*> 

_  v2a2BR 
2R?R1 

//,* 


4X  4Y 


+  4<r2y2(  1-^)1  d£  dri 


00 

f  2 

{  CXP  ~ 

r  2  2 

—7  +4VV(l-^) 

^  / 

rdr 


(7) 


(8) 
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For  general  forms  of  4*  »  Eqs.  7  and  8  are  rather  difficult  to  integrate 

and  it  has  been  customary  to  examine  the  high-  and  low-frequency  cases 
separately: 

high  frequency  »  1 

low  frequency  4y^cr^  «  1 

Perhaps  we  can  illustrate  the  difference  between  the  high  and  low 
cases,  or  the  smooth  and  rough  bottoms  by  referring  to  Fig.  6.  The 
upper  profile  was  taken  over  an  area  presumed  to  be  rough.  The  dis¬ 
similarity  of  the  reflected  signal  from  one  ping  to  the  next  is  quite  evident. 
Although  we  can  follow  reflection  interfaces,  the  reflection  signals  ap¬ 
pear  somewhat  sporadically.  The  lower  profile  was  taken  over  a  smooth 
area  and  here  the  signals  from  the  reflection  interfaces  are  much  easier 
to  follow. 


Fig.  6.  Continuous  seismic 
profiles  taken  over  smooth 
and  rough  bottoms. 
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The  reader  may  have  noticed  that  we  have  said  very  little  about 
the  reflection  coefficient.  We  assume  R.  is  the  reflection  at  the  water 
bottom  interface  and  can  include  subbottom  layers.  It  is  given  in  the 
appendix  and  is  discussed  in  more  detail  by  Tolstoy  and  Clay  (Ref.  4, 

Sections  2.  4  and  2.  5). 

5.  The  Illumination  and  Bottom  Roughness 
The  roughness  of  the  bottom  can  be  exceedingly  complex, 
cr^  and  ^  are  unknown  and  are  the  quantities  that  we  are  hoping  to 
determine.  As  we  mentioned  earlier,  our  procedure  is  to  make  a  guess, 
test  it,  and  then  improve  our  guess.  A  few  remarks  about  the  form  of  ^ 
for  several  expected  types  of  bottom  should  be  helpful.  In  this  it  is 
important  to  recognize  that  the  illumination  function  limits  our  field  to 
the  illuminated  area  and  that  we  are  considering  only  the  roughness  over 
the  illuminated  area.  As  shown  in  Fig.  7,  roughness  outside  the  illuminated 
area  does  not  matter  —  at  least  so  long  as  it  does  not  cause  shadows.  A 
bottom  that  appears  to  be  moderately  rough  on  an  echo  sounding  profile 
could  be  smooth  over  the  illuminated  area  as  sketched  in  Fig.  8.  These 
examples  show  that  we  are  concerned  about  bottom  irregularities  that  have 
wavelengths  less  than  the  dimensions  of  the  illuminated  area,  X  and  Y 
Correlation  functions,  vj'  *  are  sketched  in  Fig.  9  for  typical 
bottom  types.  The  correlation  function  can  be  regarded  as  being  zero  for 
distances  greater  than  L  .  We  shall  use  these  properties  in  our  ap¬ 
proximate  evaluations  of  the  integral. 
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Fig.  7.  The  illuminated  area.  Fig.  8.  Smooth  bottom  over  A  , 

The  area  is  indicated  by  the  the  illuminated  area, 

shading  on  the  bottom  profile. 


L  f 


TV>f>?n??rtr?frV7f>' 


Fig.  9.  Bottom  profile  and  correlation  functions 


6,  High-Frequency  Signals,  (■yc)^ 

Much  of  the  ocean  has  a  bottom  that  can  be  considered  rough  for 
most  of  the  frequency  range  used  in  geophysical  and  sonar  applications. 
Generally  speaking,  the  bottom  is  rough  when  the  rms  roughness  is 
many  acoustical  wavelengths. 

A  number  of  experiments  have  shown  that  the  sound  is  scattered 

% 

incoherently  and  we  can  then  recast  the  form  of  the  equation  for  (pp  ) 
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It  is  convenient  to  rewrite  Eq.  (8)  so  that  the  integral  is  proportional  to 
a  scattering  function.  In  doing  this,  we  are  assuming  the  validity  of  the 
incoherent  scattering  interpretation.  With  the  aid  of  Eq.  (8),  let  us  define 
the  high-frequency  scattering  function  and  other  quantities  as  follows: 


y2ffi2 

2ir 


QO 

£  D  exp[-4y2cr2(l -40]  rdr 


for 


y2<r2  »  1 


(9) 


<PP*>  =  p2AS/R2  (10) 

Pi  =  B2/R2  (11) 

A  =  irR2  (or  XY)  .  (12) 


Our  purpose  in  introducing  S  is  to  separate  the  scattering  function  so 
that  we  can  display  more  easily  its  dependence  upon  frequency,  angle,  and 
4*  .  On  referring  to  Fig.  3,  S  can  also  be  written  as 

s  ■  (PpVfth,  s)  2 

(P.h,s)2  *  »1A/R2 

3 

Eckart  suggested  that  the  integral  in  (9)  can  be  evaluated  rather  easily 
by  recognizing  that  the  principal  contribution  comes  from  the  very  small 
values  of  r  and  that  the  integral  is  nearly  zero  for  r  greater  than 
the  correlation  distance  L  .  Also,  the  illumination  function  D 
assures  us  of  convergence  for  large  r  .  It  has  been  customary  to 
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assume  an  exponential  form  for  4*  and  then  to  expand  4*  as  a  power 
3  5 

series.  ’  Depending  upon  which  terms  are  kept,  one  finds  the  scattering 
function  to  have  quite  different  dependencies  upon  frequency.  For  illustra¬ 
tion,  we  calculate  S  for  the  following  simple  types  of 


4>  =  a  0  <  r  <  L 

4*  =  0  r>L  (13) 

a  <  1 


4<  ~  e'r  /2L  ,  or  41  =  l-r2/2L2  0  <  r  <  2L2 

q  q 

4/  =  0  r2  >  2L2 

4j£  ~  e  ,  or  4^  =  1-r/L  0  :£  r  S  L 

v|i£  =  0  r  >  L 


(14) 


(15) 


These  functions  are  illustrated  in  Fig.  10.  The  double  cross  hatch 
indicates  the  region  that  makes  the  main  contribution  to  the  integral. 
The  skirts  of  the  correlation  function  are  probably  unimportant. 


Fig.  10.  41  for  high-frequency  approximation. 
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The  dependence  of  D  upon  r  has  also  been  sketched  on  ij/ 

>.o  show  that  D  is  nearly  one  over  the  range  of  values  of  R  for  which 
is  different  from  zero.  We  can  make  an  approximate  evaluation  of  (9) 
by  assuming  that  D  is  one  over  the  range  of  r  from  zero  to  L  . 
The  value  of  the  integral  over  D  for  X  =  0  gives  an  estimate  of  the 
magnitude  of  the  errors  involved  in  the  approximations  as  follows: 


y2®2 

2tt 

00  .22 
/  D  "  rdr 

0 

(16) 

y2®2 

4tr 

,  .  2  2 

R2  1 

(17) 

For  y2<r2  =  10  ,  the  exponent  in  Eq.  17  is  of  the  order  of  e~^ 

7  16 

The  other  factors  in  (17)  are  probably  less  than  10  (or  e  )  for 
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R  -  1000  X  ,  and  (17)  is  less  than  e  or  about  -100  dB.  The  scattering 
functions  are  the  following: 


4\ 


+1  ' 


S.  * 


s  * 


s.  * 


4ir 

a2  l2 
TET  ~l 


y2L2  exp  -[4y2<r2(l-a)] 


A 


„  2  2 
32tt  y  <t 


$ 


(18) 

(19) 

(20) 


The  dependence  of  these  three  functions  upon  grazing  angle  is 
shown  in  Fig.  11.  For  illustration,  the  following  parameters  were  chosen: 


y2o-2  =  10 

*b  =  0.9 

-3 

VI 

VI 

© 

o 

II 

b 

o 

II 

?- 

• 
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Fig.  11.  Comparison  of  the  high- 
frequency  scattering  functions  in 
specular  direction  as  a  function  of 
grazing  angle,  (y ^<*"2  =  10, 

L/cr  =  10,  <J>b  =  0.  9  .) 


The  scattering  function  in  the  specular  direction  is  independent  of 
angle  and  frequency  for  quadratic  type  approximations  to  4*  ,  Eqs.  (14) 

and  (19).  Other  approximations,  such  as  the  box  and  linear  functions, 
show  angle  and  frequency  dependence.  The  type  of  correlation  function 
also  changes  the  magnitude  of  the  scattering  function,  even  though  the 
values  of  <r  and  L  are  the  same. 

The  scattering  of  sound  into  directions  other  than  the  specular 
also  depends  upon  the  type  of  correlation  function.  For  non-specular 
scattering,  we  calculate  the  mean  square  of  Eq.  (1)  and  reduce  the 
multiple  integrals.  After  some  manipulation,  one  obtains  the  following 
for  the  isotropic  case: 

2  2  00 

S  “  J  Jo  (2Kr)  exp[-4Y2<r2(l-*)]  rdr  (21) 

where  D  has  been  replaced  by  1  in  view  of  the  rapid  convergence 
of  the  integrand.  The  scattering  functions  for  4^  and  4^  are: 
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(22) 


S 


t 


2.2  2  2 
y  L  y  (T 

“  4  4  2.  2.  3/2 

(4v  ir  +k  L  )  ' 


(23) 


For  the  general  case,  y  and  k  are  functions  of  the  directions 
to  the  source  and  receiver  as  shown  in  Fig.  5.  High-frequency  scattered 
sound  has  been  studied  extensively  by  C.  W.  Horton  and  others  at  the 
University  of  Texas.  The  comparison  of  theory  and  experiment  shown  in 
Fig.  12  was  taken  from  Horton  et  al.  ^  In  a  preceding  paper,  Horton  and 

5 

Muir  show  that  some  modification  of  the  surface  boundary  condition  is 

needed  to  improve  the  theoretical  fit  to  the  data.  Their  modification 

appears  to  be  approximately  a  shadowing  factor.  Their  expression  for 

S„  is  (in  our  notation) 
l 


mr_ 

1 6ir 


2.2  2  2 

y  L  y2«r 

hAVlV/2 


(23a) 


■y ^  =  k  cos 


The  correlation  function  of  their  rough  surface  had  an  approximately 
linear  form  for  small  r  .  The  example  we  have  chosen  is  for 
vertically  incident  sound. 
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r  Fig.  12.  The  scattering  coefficient 

for  the  high-frequency  case.  b  The 
dashed  line  is  Eq.  24a.  f=100  kHz  , 
=0.92  cm  ,  L=6.  5  cm  ,  k<r=3.  9  . 
The  receiver  angle  is  grazing  angle 
and  the  arrow  indicates  the  source 
position,  90*  grazing  angle. 


For  a  second  illustration,  we  use  the  so-called  reverberation  and 
determine  the  back-scattering  functions,  i.  e.,  .  The  back- 

scattering  functions  for  and  4^  are  shown  in  Fig.  13.  The  strong 

scattering  of  sound  in  the  specular  direction  (<{>  =  90*  grazing  angle)  is 
quite  evident  for  the  bottom  having  the  quadratic  correlation  coefficient. 

A  bottom  having  a  linear  correlation  coefficient  scatters  sound  much  more 
diffusely. 


Fis.  13.  Back-scattering  functions. 
(yZcr2  =  io  ,  L/(r  =10  . ) 
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The  bottom  scanning  sonar  "looks"  at  the  back  scatter  and  changes 


of  bottom  type  may  show  as  high  or  low  scattering  areas.  The  lateral 
echo  sounder  version  is  shown  in  Fig.  14.  The  transducer  transmits  a 
short  ping  and  the  travel  time,  T  ,  of  the  returning  sound  depends 
upon  range  and  angle.  The  intensity  of  the  signal  is  given  by  Eq.  (10) 
and  is  approximately  the  following: 


<p2)  ~  t'4as 

where  A  ~  TAT  (angular  beam  width) 

2  _3 

Since  (p  )  decreases  as  T  ,  it  is  evident  that  the  receiver  must 

have  a  very  wide  range  of  time-varying  gain  and/or  automatic  gain  control. 

The  instruments  are  usually  operated  about  one-fifth  their  maximum  range 

above  the  bottom  and  the  variation  of  (p2)  would  be  expected  to  be 
4 

larger  than  10  . 
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Two  examples  of  lateral  echo  sounder  records  are  shown  in  Fig.  15. 
The  featureless  record  was  taken  in  a  smooth  area  and  the  other  was  taken 


in  a  rougher  area.  The  rough  area  shows  changes  of  back  scatter  that  are 
probably  due  to  changes  of  slope.  Light  areas  or  shadows  are  usually  to 
the  left. 


7 

Fig.  15.  Lateral  echo  sounder  records. 

(Darker  is  stronger  back  scatter. ) 

In  summary,  the  dependence  of  the  high-frequency  scattering 
coefficient  upon  frequency  and  angle  can  be  used  to  guess  L/r  and  the 
dependence  of  ^  upon  r  for  small  r 

7.  Low-Frequency  Signals,  y0"  <  1 
Studies  of  acoustical  experiments  performed  over  relatively  smooth 
bottoms  have  shown  that  the  energy  is  largely  reflected  and  that  only  a 
small  fraction  is  scattered.  Since  the  Helmholtz  equation  yields  a  solution 


1 


that  is  superficially  similar  to  the  incoherent  scattering  solution,  we 
think  it  desirable  to  commence  our  discussion  of  the  reflection  interpreta¬ 
tion  by  considering  Eq.  (1). 

We  follow  a  line  of  reasoning  suggested  by  Eckart  and  discussed 

4 

in  detail  by  Tolstoy  and  Clay.  Let  us  assume  that  an  acoustic  interaction 
experiment  such  as  shown  in  Fig.  1  is  repeated  many  times.  But  in 
between  each  repetition,  the  black  box  is  changed  so  that  each  received 
signal  is  different.  If  the  properties  of  the  black  box  are  assumed  to 
have  a  statistical  distribution,  then  one  can  average  over  the  distribution 
function. 

Let  us  assume  that  £  ,  the  roughness,  has  a  Gaussian  distribution 
over  the  illuminated  area.  The  average  of  p  over  a  surface  having  an 
rms  roughness  <r  is 

ik(R  +R  )  ,  2 

<p>  ■  -iB^felRz  e'2v  '  •  (24) 

Notice  that  for  <r  =  0  ,  or  a  perfectly  smooth  surface,  the  form  of  (24) 

remains  the  same.  Equation  (1)  or  (24)  for  £  =  0  or  c  =  0  is  the 
Helmholtz  way  of  expressing  the  signal  reflected  when  the  surface  is 
illuminated  by  means  of  a  directional  source  (6a). 

P|  Q  =  -gBXy  R  exp  {-2u2XZ-2(32Y2fik(R1fR2)}  .  (24a) 

: 

|  Following  the  idea  shown  in  Fig.  2,  it  is  convenient  to  define  the  reflection 

signal  reflected  by  a  perfect  reflector  as 
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pth.  r  =  <p> 

tn’ r  <r=o,  m=i 


and  then  express  (25)  as  follows: 


<P>  =  Pth,r<»> 


(a)  =  a 


-.2  2 
-2\  <r 


The  coherent  reflection  coefficient  describes  the  reflection  in  the  specular 
direction.  Since  the  signal  from  an  omnidirectional  source  can  be  ex¬ 
panded  in  terms  of  outgoing  plane  waves,  it  is  evident  that  the  coherent 
reflection  coefficient  can  be  applied  to  the  omnidirectional  source.  We 
can  follow  this  line  of  reasoning  and  assert  that  p^  r  can  be  calculated 
by  the  usual  image  method  as  shown  in  Fig.  2. 

The  coherent  reflection  coefficient  is  independent  of  the  surface 
correlation  function.  The  mean  square  signal  consists  of  two  terms,  the 
reflected  sound  and  the  scattered  sound.  The  latter  is  dependent  upon  the 
surface  correlation  coefficient  and,  as  we  shall  show  for  a  simple  case, 
measurements  of  scattered  sound  can  be  used  to  determine  the  spatial 
spectrum  of  the  surface. 

For  low  frequencies  and  small  y<r  ,  the  dependence  upon  y<r 
and  ^  can  be  approximated  as  follows: 


exp  - [4y^<r^(l -^)]  —  (l  +  4y^o’^)  e  ^ 
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We  substitute  (28)  into  (7)  and  (8)  and  find  that  the  first  term,  or  the 
multiple  integral,  is  independent  of  4*  and  is  the  square  of  the  mean 
signal,  (24)  or  (26).  The  second  term  involving  ip  is  the  scattered 
sound.  On  making  these  substitutions  and  recalling  (10)  -  (12),  we  obtain 


<PP*>  -  <P>2  + 


2  2m2 
V  O’  JR 
“2“ 
ttR, 


P2  A  /  / D"iP(e,  T!)  e2i(a^+PTl)dedr1  .(29) 


If  ip  vanishes  at  values  of  £,  q  much  smaller  than  X  and  Y  ,  then 
again  the  influence  of  D  on  the  integral  can  be  ignored.  The  integral  is 

5 

then  the  two-dimensional  Fourier  transform  of  the  space  correlation  function. 

It  can  be  written  as  follows: 

*  =  7‘ P(|,q)  e2i(a|+pTl)dedT1  (30) 

-oo 

where  is  the  power  spectrum  of  the  surface  roughness.  The  substitu¬ 

tion  of  (30)  into  (29)  yields 

(pp")  -  <P>2  +  \2(r2  +  •  (31) 

ttR, 


The  second  term  represents  the  incoherent  scattered  sound  and  has  the 
same  geometrical  dependencies  as  given  in  Fig.  3.  Equation  (31)  is  also 
the  acoustical  analogue  of  the  familiar  optical  problem  of  the  diffraction 
grating.  Here  we  use  a  known  continuous  wave  signal  to  determine  the 
spectrum  of  the  rough  surface  rather  than  a  known  grating  to  analyze  the 
spectrum  of  an  unknown  signal,  vp  and  the  corresponding  V  for 
several  simple  functions  are  the  following: 
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+  1  =  e-<e2+n2)/2L2 


lirLZ'-2{*W)LZ 


+2  .  e-'«2t”2»/2L2cosa6  *, 

+  3  =  e  *3 

for  L  «  X  or  Y 


ttLi2{ 


2;e-(a±2a)2L2/2  } 


r-  2  -Z^L\ 

2\f2v  L  e _ 

l+4a2L2 


The  correlation  functions  and  corresponding  spectrum  functions 
are  sketched  in  Fig.  16.  Numerical  examples  of  the  total  and  scattered  sound 
field  as  a  function  of  a  are  shown  in  Fig.  17.  The  width  of  the  reflected 
sound  field  is  dependent  upon  the  dimensions  of  the  illuminated  area.  That 
is,  the  larger  the  illuminated  area,  the  narrower  the  reflected  sound  beam. 
The  width  of  the  scattered  sound  field  varies  similarly  with  the  correlation 
distance. 


Fig.  16.  Correlation  functions 
and  power  spectrum  functions. 


X  =  2  0  m 


Fig.  17.  Scattered  and  reflected 
sound.  The  surface  spectrum  is  ^  . 
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In  Fig.  18  we  show  the  scattered  sound  for  a  peaked  spectrum 
function  ^  as  sketched  in  Fig.  16.  The  specular  reflection  is  very 
much  narrower  than  the  scattered  sound  because  is  10  .  The 

width  of  the  scattered  sound  peaks  would  be  about  the  same  size  as  the 
reflected  sound  for  . 

S"«  #* 


Fig.  18.  Reflected  and  scattered  Fig.  19.  Theoretical  and  experi- 

sound  for  V  .  mental  scattering  of  sound  for  small 

2  2  8 
Y  a 

The  sound  scattered  by  a  surface  having  a  spectrum  of  the  type  flr 

was  studied  by  means  of  laboratory  experiments  at  Brown  University. 

8 

The  theoretical  curve  and  data  of  Proud,  Beyer,  and  Tamarkin  are 
shown  in  Fig.  19. 

We  should  say  a  few  words  about  experiments  to  measure  ^  . 

First,  the  scattered  sound  is  proportional  to  (y^cr^)  and  that  was 
assumed  to  be  small.  Thus,  the  ratio  of  the  scattered  sound  to  the  co¬ 
herently  reflected  sound  is  going  to  be  very  small.  If  we  were  to  use  an 
omnidirectional  source,  then  whatever  the  direction  in  which  we  make  our 
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observation,  we  will  predominantly  receive  the  coherently  reflected 
component.  If  we  confine  our  experiments  to  single  source  and  single 
receiver  experiments,  then  the  best  that  we  can  do  is  to  measure  <r  . 
Perhaps  we  can  separate  the  coherent  component  from  the  scattered 
component  by  the  use  of  receiving  arrays  and  thus  measure  ^  . 
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PART  II.  INTERPRETATION  OF  MARINE  GEOPHYSICAL  DATA 


1.  Introduction 

Preliminary  examination  of  the  marine  geophysical  data  indicates 
that  the  magnitude  of  the  bottom  loss  measurements  can  be  related  to 
physiographic  provinces  and  types  of  ocean  bottom.  As  we  shall  show 
later,  these  measurements  can  be  used  to  estimate  some  of  the  important 
parameters  that  describe  the  bottom  of  the  ocean.  The  procedure  that  we 
suggest  for  the  interpretation  of  reflection  data  combines  our  knowledge 
of  the  unconsolidated  marine  sediments  and  reflection  of  signals  at  a  rough 
interface.  It  will  be  shown  that  roughness  of  the  bottom  frequently  plays  a 
dominant  role  in  the  areas  where  very  little  of  the  acoustical  energy  is 
reflected  at  the  bottom  of  the  ocean. 

A  physical  approach  is  needed  because  a  bottom  loss  measurement 
does  not  have  general  meaning  or  applicability  in  itself.  The  loss  measure¬ 
ments  must  be  considered  in  the  context  of  the  experiment  and  other 
geological  and  geophysical  data.  (Although  we  shall  not  pursue  the  matter, 
a  loss  measurement  may  be  applicable  only  to  equipment  that  works  exactly 
in  the  same  way  as  the  measurements  were  made.)  Our  purpose  in  this  is 
to  supply  some  of  the  context  and  relationships  that  can  give  meaning  to 
loss  measurements.  Perhaps  it  is  needless  to  say  it,  but  knowledge  of 
the  ocean,  its  bottom,  and  acoustical  theory  can  be  combined  to  predict 
sonar  performance. 


2.  Summary 

Much  of  the  deep  ocean  is  covered  by  sediments  that  can  be  described 
as  sands  and  clays.  Often  the  layers  are  a  fraction  of  a  meter  thick  and 
alternate.  The  range  of  the  values  of  the  average  velocities  and  densities 
of  the  top  few  meters  of  sediments  is  such  that  the  maximum  range  of  re¬ 
flection  loss  at  normal  incidence  is  approximately  7  to  14  dB.  If  the  bottom 
is  known  to  be  sand  or  denser  sediment,  the  losses  would  be  expected  to 
range  from  7  to  10  dB.  Inasmuch  as  low-frequency  acoustic  signals  have 
longer  wavelengths,  these  signals  are  reflected  by  the  deeper  more  com¬ 
pacted  layers  and  the  reflection  losses  are  smaller.  The  reflection  from 
a  layered  smooth  bottom  is  both  frequency  and  angle  dependent.  A  brief 
discussion  is  given  in  Appendix  A.  It  would  appear  that  bottom  losses 
larger  than  10  dB  and  certainly  those  larger  than  14  dB  must  be  due  to 
some  mechanism  other  than  transmission  into  the  bottom. 

The  results  of  the  marine  geophysical  surveys  indicate  that  the 
reflection  losses  are  high  in  areas  where  the  bottom  is  known  to  be  rough 
and  the  losses  are  small  in  smooth  areas.  The  physiographic  charts  such 
as  those  of  Heezen  and  Tharp^are  particularly  useful  in  determining  the 
location  of  the  smooth  and  rough  areas.  To  be  useful  for  sonar  perform¬ 
ance  prediction,  the  description  of  the  bottom  needs  to  be  much  more 
precise  than  smooth  or  rough.  In  the  following,  we  give  guide  lines  for 
quantitative  interpretation  of  the  marine  geophysical  reflection  measure¬ 
ments. 
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a)  Smooth  and  Slightly  Rough  Ocean  Bottom 
Over  a  smooth  bottom  at  vertical  incidence,  8  to  10  dB  losses  were 
measured  at  3.  5  kHz  and  slightly  smaller  dB  losses  were  measured  at 
0.  5  kHz.  Even  in  the  smoothest  areas  we  frequently  find  that  the  vertical 
incidence  bottom  loss  at  3.  5  kHz  is  a  little  larger  than  that  for  0.  5  kHz. 
Also  in  smooth  areas,  the  seismic  profiles  show  smooth  layered  initial 
bottom  and  subbottom  reflection  as  shown  in  Fig.  20.  The  precision  echo 
soundings  (12.  5  kHz)  usually  show  good  bottom  reflection  and  some  sub¬ 
bottom  reflections. 


o 

x> 


Fig.  20a.  Smooth  bottom, 
<r  <  0. 1  m  (r  is  the  rms 
roughness). 


Fig.  20b.  Sketch  of  seismic 
profile. 


As  shown  in  Fig.  20a,  the  dependence  of  bottom  loss  on  grazing 
angles  is  about  that  expected.  In  smooth  areas,  0.5  and  3.5  Hz  reflection 
coefficients  of  the  ocean  bottom  tend  to  0  dB  at  small  grazing  angles.  This 
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is  a  typical  result  and  is  based  upon  the  Marine  Geophysical  Survey  data. 

The  method  used  to  estimate  cr  is  discussed  later  and  in  Appendix  B.  As 
part  of  the  routine  geophysical  surveys,  bottom  losses  are  measured  with  a 
3.  5  kHz  echo  sounder.  These  data  show  about  8  to  10  dB  loss  in  smooth 
areas  and  appear  to  compare  reasonably  well  with  the  acoustic  reflectivity 
surveys.  We  mentioned  earlier  that  the  reflectivity  correlates  with  the 
physiographic  province.  Measurements  indicate  that  abyssal  plains  vary 
from  very  smooth,  <r  <  0.  1  m  ,  to  slightly  rough,  <r  <  0.  5  m  . 
b)  Intermediate  Roughness  of  Bottom 

Areas  having  an  intermediate  roughness  have  bottom  losses  that  are 
frequency  dependent.  As  shown  in  Fig.  21,  the  loss  at  3.  5  kHz  is  considerably 
higher  than  the  loss  at  0.  5  kHz.  The  seismic  profiles  are  essentially  a  graphic 
display  of  vertical  incidence  reflection  data.  The  seismic  signals  are  generally 
low-pass  filtered  to  pass  less  than  0.  2  kHz  and  as  such  are  an  indication  of  the 
0.2  kHz  reflection  coefficient.  In  these  areas,  the  seismic  profiler  data  are 
good,  i.  e.,  the  bottom  and  subbottom  reflections  are  well  defined,  but  the 
subbottom  structure  changes  over  short  distances. 

The  amount  of  scattering  depends  upon  the  ratio  of  the  roughness  to 
the  acoustic  wavelength.  In  areas  having  tr  -  1  m  ,  the  3.  5  kHz  signal  has 
such  a  short  wavelength  that  the  sound  is  scattered  at  all  grazing  angles. 

Since  the  acoustic  wavelength  of  0.5  Hz  is  seven  times  larger,  the  same 
scale  of  roughness  does  not  scatter  the  0.  5  Hz  signals  nearly  as  much.  The 
3.  5  kHz  vertical  incidence  losses  range  from  14  to  20  dB. 

Some  basins  and  many  continental  rises  have  intermediate  roughness, 

<r  =*  1  m  . 


-29- 


CD 


20 


yywwvM 


'3.5  kHi 


30* 


60* 


90*^ 


Fig.  21a.  Intermediate  roughness 
o-  =«  1  m. 


Fig.  21b.  Seismic  profile. 
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c)  Rough  Areas 

Rough  areas  appear  to  scatter  the  0.  5  and  3.  5  kHz  signals  about  the 
same  amount,  thus  the  bottom  losses  have  the  same  range.  Typical  reflec¬ 
tion  losses  are  sketched  in  Fig.  22a.  Several  types  of  seismic  profiles 
that  are  obtained  in  rougher  areas  are  sketched  in  Fig.  22b.  The  shaded 
portion  indicates  the  "basement"  reflection.  The  scattering  of  sound  in 
very  rough  areas  such  as  ridges  is  primarily  dependent  upon  the  mean 
square  slopes  of  the  irregularities  and  independent  of  frequency.  The 
fluctuations  would  be  expected  to  be  and  are  very  large  for  small  changes 
of  geometry. 
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Fig.  22a.  Rough  area,  cr  >  2m. 


Fig.  22b.  Typical  seismic  profiles  in  rough  areas. 
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3.  Formulas  and  Discussion  of  Theory 

The  reflection  and  scattering  theory  is  given  in  detail  by  Tolstoy 

and  Clay  (Ref.  4,  Chapter  6).  Results  are  summarized  in  the  Appendices. 

In  this  section  we  give  a  few  results  and  rules  of  thumb  to  estimate  bottom 

roughness  from  bottom  loss  measurements.  Some  aspects  of  the  measure- 

9 

ment  of  roughness  are  discussed  by  Clay. 

a)  Bottom  Reflection  -  Smooth  Area 
As  shown  in  Fig.  2,  Part  I,  the  reflected  signal  p  ^  is»  with 
the  aid  of  the  Eckart  reflection  coefficient^  (Appendix  Eqs.(Bl)and(B2)) 
the  following.  (The  reflection  coefficient  R^  or  ®-  is  given  by  Eqs. 
(Al),  (A5),  ,  or  (All)  in  the  Appendix.  The  choice  of  expression  depends 
upon  how  complicated  the  bottom  is.) 


where 


2  Po<®>2 

Pr6f  '  <V*2>2 
,  2_2 

(JB.)  =  fflL  e"2v 


(35) 


(36) 


p  is  acoustical  pressure  at  unit  distance 
y  =  ^  sin  a  (37) 

y  is  the  vertical  component  of  wave  number 


The  exponential  factor  is  e  for 

a-  =  —1 -  (38) 

trsIZ  sin  4> 
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and  (38)  can  be  used  as  a  criterion  of  roughness.  The  bottom  loss  in  dB 
can  be  written  with  the  aid  of  Eq.  (Dl)and  is 


-20  log10(S>  =  -20  log  Q  JR  +  680  ^  sin2  <>  .  (39) 

\ 

Equation  (39) can  be  used  to  estimate  the  roughness  of  the  area  as  follows: 

At  signal  grazing  angles  JR.  tends  to  1,  thus  the  bottom  loss  is  ap¬ 
proximately  the  following: 

tr2  2 

[Bottom  loss  for  small  4>]  —  680  — j  sin4^  .  (40) 

The  value  of  <r  can  be  calculated  with  the  aid  of  Eq.  (40)  for  reflection 
coefficients  that  depend  upon  <J>  as  sketched  in  Figs.  20  and  21. 

For  a  given  roughness,  it  is  evident  that  the  reflection  loss, 

Eq.  (37),  is  the  same  for  a  100-Hz  signal  at  30°  grazing  angle  and  a 
500-Hz  signal  at  6".  That  is 

at  100  Hz  y  =  —  100  sin  30°  =  —(0.5)  (100) 

c  c 

(41) 

at  500  Hz  y  =  —  500  sin  6#  =  —  (0.  1)  (500) 

c  c 

b)  Bottom  Scattering  -  Rough  Area 
In  rough  areas  as  shown  in  Fig,  22,  most  of  the  energy  is  scattered. 

3 

Under  these  conditions,  as  shown  by  Eckart  and  in  Appendix  C,  the 
scattering  is  dependent  upon  the  mean  square  slopes  or  <r  and  the 
correlation  distance.  The  scattered  sound  is 
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(42) 


sc 


-  m.  sa 


R1R2 


A  is  the  scattering  area 
S  is  the  scattering  function 


Comparison  of  (35)  and  (42)  shows  that  the  reflected  and  scattered 
sounds  have  different  functional  dependence  upon  range.  The  reflected 

signal  depends  upon  range  as  (R^+R^)*'2  ,  whereas  the  scattered  signal 

2  _  2  -  ? 

p  depends  upon  range  as  R.  R  “  .  Calculation  of  the  scattered 

SC  1  b 

sound  requires  knowledge  of  the  scattering  area  A  .  Note  that  in 
many  experimental  measurements  A  has  not  been  very  well  defined 
and  only  an  apparent  reflection  coefficient  has  been  determined.  The 
detailed  discussion  of  a  means  of  estimating  A  for  particular  experi¬ 
ments  is  given  in  Appendix  C. 


4.  Conclusion 

The  marine  geophysical  data  should  be  used  with  great  care.  In 
smooth  areas,  the  experimental  data  can  be  used  to  estimate  the  bottom 
parameters  and  the  roughness.  Our  ability  to  use  surface  reflection  to 
make  meaningful  estimates  of  the  sound  velocity  in  the  bottom  layers 
decreases  as  the  roughness  increases.  In  rough  areas  we  must  know 
the  scattering  area  in  order  to  estimate  the  scattering  coefficient. 
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APPENDIX  A 


REFLECTION  AT  A  SMOOTH  BOTTOM 

Much  of  the  ocean  bottom  is  covered  by  sediments  that  can  be  de¬ 
scribed  as  clays  and  sands  and  often  these  layers  alternate.  The  acoustic  ve¬ 
locity  in  the  sediments  shows  an  average  increase  with  subbottom  depth.  The 
temperature  gradient  in  the  subbottom  also  causes  an  increase  of  acoustic  ve¬ 
locity.  Nafe  and  Drake*^  estimated  the  gradient  to  be  from  0.4  sec-*  to  1.0 
sec  *.  Depending  upon  the  material,  the  acoustic  velocity  at  the  water  sedi¬ 
ment  interface  may  vary  from  a  few  percent  less  than  that  in  the  water  to  some¬ 
what  higher.  It  is  well  known  that  the  acoustic  velocity  in  fluid-saturated  uncon¬ 
solidated  sediments  also  depends  upon  the  temperature.  Because  of  the  thermal 
gradient  in  the  sediments,  the  acoustic  velocity  would  also  be  expected  to  have 
an  average  increase  as  a  function  of  depth.  Laboratory  acoustical  measure¬ 
ments  of  recovered  cores  should  be  made  at  the  proper  temperatures.  Data 
taken  at  some  locations  in  the  Pacific  Ocean  indicate  that  the  acoustic  velocity 
in  the  sediment  near  the  water-sediment  interface  is  less  than  the  acoustic  ve¬ 
locity  in  the  water.  Measurements  of  the  acoustic  velocity  of  the  sediments  in 
the  Hatteras  abyssal  plain  in  the  Atlantic  Ocean  have  been  around  1.6,  1.7 
km/sec.  **  On  the  average,  the  acoustic  velocity  in  the  sediments  increases 
with  subbottom  depth.  The  sound  ray  paths  for  both  cases  are  shown  in  Fig.  Al. 

Experimental  measurements  of  the  attenuation  of  acoustic  waves  in 
sediments  have  shown  that  the  materials  absorb  energy.  The  absorption  is 
dependent  upon  the  kind  of  sediment  and  frequency.  Generally  the  absorption 
is  larger  at  higher  frequencies.  In  addition  to  roughness,  both  the  layering 
and  absorption  cause  the  reflection  coefficient  to  be  dependent  upon  frequency. 

In  view  of  this,  unqualified  attempts  to  relate  experimental  data  to  a  simple 


non-layercd  model  may  be  misleading.  The  conversion  of  compressional 
wave  energy  to  shear  waves  can  also  affect  the  reflection  coefficient. 


WATER 

/ / //7//Z // /  SEDIMENT 

<p  <  CRITICAL 


Fig.  Al.  Comparison  of  ray  paths  for  higher  and 
lower  velocities  of  sound  at  the  water-sediment  interface. 

We  assume  that  the  reader  is  familiar  with  the  reflection  from  the 
plane  interface  between  mediums  1  and  2  and  the  conditions  for  vertical 
reflection.  For  the  geometry  shown  in  Fig.  A2,  the  reflection  coefficient 
Rj2  as  given  by  Tolstoy  and  Clay  (Ref.  4,  p.  30)  is 


P2Vl"PlV2 

P2VP1^2 


(Al) 
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where 


v,  =  fj-sin*!  -  kjSin^, 

y2  =  f-»ln+2  =  k2Sin*2 


and  Snell' s  Law  gives  the  relation  between  and 


cos  4>j  cos  4>2 


q  =  angular  frequency 
c ^  =  velocity  in  medium  1 

c y  =  velocity  in  medium  2 

M 

and  are  the  densities. 


(A  2) 


(A  3) 


Fig.  A2.  Simple  reflection. 
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Critical  reflection  occurs  for  grazing  angles  less  than  <£ 


lc 


lR  12 1  =  1  for  $1  <  ‘f'lc 


cos  ^ lc 


(A4) 


The  usual  form  of  the  reflection  equation  can  be  obtained  by  the  substitution 
of  (A2)  into  (Al)  and  the  elimination  of  w  .  The  reflection  from  a  layer 
is  frequency  dependent  and,  from  Tolstoy  and  Clay  (Ref.  4,  p.  30),  it  is 


1R 


2iy  h_ 

R12  *  R23  6 _ 

2iV2h2 

l  +  R.  -R.,  e 


h^  is  layer 
c  thickness. 


(A  5) 


The  reflection  coefficient  for  a  thin  layer  is  shown  in  Fig.  A3.  The 

maxima  and  minima  are  caused  by  the  interference  of  the  reflections  from 

the  two  interfaces.  Although  roughness  is  treated  later,  the  dashed  curve 

shows  the  effect  of  roughness  on  the  coherent  reflection  coefficient. 

Multiple  subbottom  layers  can  be  treated  by  the  recursive  technique 

4 

given  in  Tolstoy  and  Clay. 
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1.5  km/sec 


Fig.  A3.  Reflection  from  a  layer.  The  solid  curve  is  the  reflection 
from  a  uniform  layer  =  <r2  =  ^  dashed  curve  is  the 

reflection  from  a  slightly  irregular  layer  <r,  =  <r,  =  1  m 


For  grazing  angles  less  than  the  critical  angle  at  the  2-3  interface, 


is  complex  and  R^  can  be  written  as  follows 


R23  =  6 


iX 


a->  <  ao 

2  2  cnt. 


(A6) 


The  substitution  of  (A6)  into  (A5)  and  calculation  of  the  absolute  square  of 
(A5)  yields  an  extremely  important  result 

|m|  =  1  .  (A  7) 

Thus  the  effective  reflection  coefficient  of  the  layer  is  one,  even  though 
the  top  layer  may  have  an  acoustic  velocity  less  than  that  in  the  water 
above  the  interface.  In  view  of  the  observation  that  the  sound  velocity  in 
sediments  increases  with  depth  of  burial,  one  would  generally  expect  to 
observe  total  reflection  phenomena  at  shallow  grazing  angles. 

Studies  of  reflection  and  transmission  in  deep  sea  sediments  have 
shown  that  the  absorption  loss  of  acoustical  energy  in  the  sediments  is 
considerably  larger  than  the  loss  in  water.  As  an  example,  we  examine 
the  characteristics  of  a  layer  having  a  small  absorption.  Following 
Tolstoy  and  Clay  (Ref.  4,  p.  221),  we  let 

\2  =  Y'2+iV2  (A8) 
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where 


y1  is  the  vertical  component  of  wave  number  in  medium  2 

w 


The  absorption  coefficient  -y"  is  roughly  proportional  to  frequency.  At 
very  high  frequency,  the  signal  reflected  from  the  second  interface  is 
highly  attenuated  and  JR  is  approximately  •  (The  reflection  from 

the  first  interface  is  also  somewhat  attenuated  by  the  lossy  bottom.) 


IR  »  Ru  (A10) 

for  ^2^2  >:>  1 

Equation  (A9)  tends  to  (A5)  at  low  frequencies,  and  we  would  expect  to 

observe  the  constructive  and  destructive  interference  of  reflections  from 

the  two  interfaces.  A  nice  set  of  laboratory  measurements  has  been 

13  14 

made  by  Barnard  et  al.  Menotti  et  al.  have  made  field  studies. 

To  conclude  this  section,  we  use  Eq.  (Al)  to  get  an  estimate  of  the 
gross  range  of  the  numerical  values  of  R^2  at  vertical  incidence. 
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At  vertical  incidence,  a  =  a  =  90*  and  R  i- 

i  c.  12' 


vert 


P 2C2-P i c j 

=  P2C2+P1C1  at  vertical  incidence*  (All) 


Since  the  c2  ranges  from  a  few  percent  below  C]  to  somewhat  above 
in  many  marine  sediments,  the  value  of  mainly  depends  upon  the 

densities.  For  a  density  range  of  1. 5  to  2.  5,  the  reflection  coefficient  is 
in  the  following  range 


0,2  <  R12  <  0,4  at  vertical  incidence 


(A  12) 


or  the  bottom  loss  ranges  from  14  to  7  dB.  In  view  of  the  average  increase 
of  pc  with  increasing  depth,  one  would  expect  to  measure  larger  re¬ 
flection  coefficients  at  low  frequencies. 
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APPENDIX  B 


REFLECTION  AT  A  SLIGHTLY  ROUGH  BOTTOM 

Irregularities  of  the  interface  cause  scattering  of  the  acoustic 
energy.  Following  the  analysis  of  Eckart,  theoretical  calculations  are 
usually  made  for  irregularity  heights  that  are  either  small  or  very  large 
compared  to  the  acoustical  wavelength  (Ref.  4,  pp.  193-220). 

For  imall  heights  of  irregularities,  much  of  the  signal  is  reflected 
in  the  specular  direction  and  the  reflection  coefficient  is  referred  to  as 
the  coherent  reflection  coefficient.  On  assuming  that  the  irregularities 
have  a  Ga’issiun  distribution  and  a  root-mean-square  roughness  <r  , 
the  reflection  coefficient  for  the  reflection  of  signals  from  a  simple  rough 
inter  {.v.  s 

,  2  2 

<*>  =  R12e '  V 

or 

(E>  =  R12e-2k2,r2sin2,|>  (Bl) 

k  =  w/c 


The  magnitude  of  the  acoustical  pressure  squared  (ignoring  the 
time  dependence)  is  given  by  the  image  method  and  is 


2 

^ref 


<VR2>‘ 


(BZ) 


where  the  geometry  is  shown  in  Fig.  Bl  and  p  is  the  pressure  at  unit 
distance.  Experimental  data  are  shown  and  compared  to  the  -ry  in  Fig.  (B2) 
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where 


<pp'  > 

2 

Pref 

♦  p2 

rsc 

<P0Po> 

2 

=  Po 

for 

tsj 

II 

2 

Pref 

<**>2 


=  1 


2 

p  is  the  sound  scattered  by  the  irregularities  on  the  surface.  This 
s  c 

energy  is  incoherent  and  while  its  intensity  is  larger  in  the  specular 
direction,  the  sound  beam  is  considerably  broadened.  The  scattered  sound 
is  usually  referred  to  as  reverberation. 


A  comparison  of  Eq.  (Bl)  and  experiment  is  the  dashed  line  in 
Fig.  B2.  The  dashed  line  is  proportional  to  (R)^  and  the  solid  line  is 
proportional  to  the  sum  of  the  coherent  reflected  signal  and  the  incoherent 
scattered  energy.  Notice  that  the  scattered  contribution  is  about  0.  1  at 
high  frequency  or  large  y <r 


-44- 


X  0*  Inc.dfnt  ingle 


-45- 


APPENDIX  C 


REFLECTION  AT  A  VERY  ROUGH  BOTTOM 


At  large  roughness  compared  to  the  acoustical  wavelengths,  the 

geometrical  dependence  of  the  scattered  radiation  is  different  from  the 

specular  case.  As  shown  in  Fig.  3,  Part  I,  the  area  A  is  assumed  to  be 

illuminated  and  the  energy  is  incoherently  scattered  by  each  part  of  the 

surface.  The  scattered  sound  p^  is 

sc 


(Cl) 


where  S  is  the  scattering  coefficient.  For  high  frequencies  and  a  very 
rough  surface  Eckart  found  that  S  is  dependent  upon  the  mean  square 
slope  of  the  surface  )  and  independent  of  frequency.  Comparison 

of  (B2)  and  (Cl)  shows  that  the  reflected  and  scattered  sound  have  different 
functional  dependence  upon  range  (Ref.  4,  p.  213).  Different  choices  of 
boundary  conditions  and  correlation  functions  of  the  irregular  surface 
lead  to  different  values  of  S  .  ^  The  effect  of  different  correlation 
functions  is  discussed  in  some  detail  in  Part  I.  Following  Eckart, 

Tolstoy  and  Clay  used  a  boundary  condition  that  assumes  the  differential 
area  to  be  a  plane  surface  and  a  Gaussian  correlation  function.  They 
found  S  to  be  the  following: 

®12 

s  at  - (Cl) 

8  "(O 

where  {tj  is  mean  square  slope. 
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Recently  Horton  et  al.  measured  the  scattered  sound  for  this 
case  and  compared  the  experimental  results  with  their  theoretical  calcula¬ 
tions.  The  correlation  function  of  their  surface  is  approximately  our 
i|;  ,  shown  in  Fig.  10,  Part  I.  A  figure  from  their  paper  is  our  Fig.  C2. 

For  these  experiments,  (iR^j  =  1  and  a  directional  source  was  used 
to  illuminate  a  known  scattering  area.  The  equation  of  the  scattering 
function  is  our  Eq.  (23a),  Part  I. 

On  the  average,  scattered  waves  add  randomly.  However,  for  a 
particular  source,  receiver,  and  rough  surface,  repeated  transmissions 
are  reproducible.  Translation  of  the  rough  surface  causes  the  scattered 
waves  to  add  differently.  A  small  change  of  the  source  or  receiver 
position  also  causes  the  waves  to  add  differently.  Thus  one  would  expect 
the  scattered  signal  level  to  have  large  fluctuations  and  the  20  dB  fluctua¬ 
tions  found  by  Horton  et  al.  would  be  expected.  The  curves  shown  in 
Fig.  Cl  are  the  average  of  five  runs.  The  fluctuations  are  present  but 
reduced. 


30  40  40  10  100  130  140 

«ECE!VE«  ANGIE 


Fig.  Cl.  The  scattering  coefficient  for  ^  =  40°  and  f  =  100  kHz. 
This  is  an  average  of  five  measurements  at  different  positions  on  the 
model  surface.  The  solid  curve  is  experimental  and  the  dashed 
curve  is  from  Eq.  8  (in  Horton  et-al.^). 
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The  data  in  Fig.  Cl  show*another  important  feature  of  the  scattered 
field.  The  scattered  radiation  is  rather  diffuse  and  slightly  peaked  in  the 
specular  direction.  The  scattering  coefficient  is  about  -20  dB.  This 
type  of  scattering  coefficient  is  used  in  Eq.  (Cl). 

If  the  bottom  is  illuminated  by  means  of  an  omnidirectional  source 
such  as  a  shot  and  the  scattered  radiation  is  received  by  hydrophone,  it 
may  be  difficult  to  determine  what  the  illuminated  area  is.  Furthermore, 
sound  reflected  or  scattered  by  subbottom  layers  may  be  present.  We 
believe  that  much  care  must  be  used  in  the  interpretation  of  data  taken  in 
rough  areas.  The  problem  is  somewhat  ill  defined.  In  addition  to 
scattering  loss  numbers,  the  readers  of  reports  should  be  given  samples  of 
the  echo  soundings,  seismic  profiles, and  graphic  records  of  the  reflection 
signals  for  each  reflection  study.  Presumably  the  signals  and  sources  of 
the  quoted  loss  numbers  would  be  identified  on  the  records. 

Let  us  return  to  equation  (Cl)  and  use  it  to  estimate  the  scattered 
sound  for  a  typical  experiment.  We  assume  that  the  source  is  a  shot  and 
the  receiver  filters,  squares,  and  averages  p^  for  At  sec.  The 
effective  illumination  area  is  essentially  given  by  the  first  At  sec 
after  the  initial  scattered  signal  arrives  (  At  »  the  source  signal  dura¬ 
tion). 

The  geometry  is  shown  in  Fig.  C2.  (The  water  is  assumed  to  be 
iso-velocity. ) 
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Fig. 

C2.  Scattering  area. 

To  determine  A  ,  we  first  estimate  and 

The  travel  times  in  the  specular  direction,  R^  =  R^  ,  are 

in  Fig.  C2. 

*1*2 

=  2R1/c1 

*1  +t2 

=  2Rj/ c 

(C3) 

At  = 

(t'l+t^)  -  (t1+t2) 

by  the  Sagitta  formula 

Y 1  “ 

^2(Atc)Rj 

(C4) 

-49- 


For  the  grazing  angle 


X1  “  VsinQ 
“  y,Vh 


The  area  of  an  ellipse  is 


A  = 


2 

1 


(C5) 


(C6) 


The  scattered  sound,  in  specular  direction  and  averaged  over  At  ,  is 
as  follows: 


sc 


=  p2(l)Hl2S-4 
At  R1 


“  p2(1)  a2  ~ 


(C  7) 


T 


~Ti 


For  constant  At  and  water  depth,  h  ,  p“  I  has  a  dependence  of 

-2  SC'At 
Rj  .  This  is  superficially  similar  to  the  ordinary  reflection  equation. 

Many  times,  reflection  coefficients  are  computed  on  the  basis  of 

the  specular  reflection  equation.  Although  it  is  incorrect,  we  proceed 


and  assume  an  apparent  reflection  coefficient 


ap 


With  the  aid  of 


Eq.  (B2)  the  corresponding  scattered  sound  in  the  specular  direction  is 


(C8) 
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The  following  value  of  IB.2 

6  ap 


is  obtained  by  equating  (C7)  and  (C8): 


“ap  =  “2s  •  <C9) 

In  the  deep  ocean  h  is  about  20,  000  ft.  On  assuming  At  —  0.01  sec, 
the  third  factor  in  Eq.  (C9)  is 

8TtAt/h  =*  Zir  X  10" 2 


and  the  apparent  reflection  coefficient  is 


B2  =*  2tr  ffi.2  S  X  10"2 
ap 


(CIO) 


As  an  example  for  scattering  in  a  very  rough  area,  let  us  assume 
that  the  rms  slope  of  all  of  the  nooks,  crannies,  and  rocks  is  15", 


<£'2)  =*  tan2 15  ‘ 


Substitution  of  the  value  into  Eq.  (C2)  yields 


S  =*  10  /  (56tt) 


(C  11) 


The  apparent  reflection  coefficient  is 


R 


2  * 
"ap 


i®2 

27 


(C  12) 


and  conversion  to  dB  yields 
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10  log  -  -14  dB  +  20  log  BSL 


(C  13) 


Near  vertical  incidence  20  leg  R  i.  of  ,he  order  of  - 10  dB  and  decreases 
to  zero  at  grazing  incidence.  Thus  we  would  expect  the  following 


10  log  ®  *  -24  dB 

ap 

-  - 14  dB 


vertical 


grazing 


in  a  rough  area. 
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APPENDIX  D 


SUMMARY 

The  dependence  of  the  scattered  sound  upon  frequency  and  the 
incident  and  reflected  angles  is  very  useful  in  deciding  whether  the 
roughness  is  greater  or  less  than  the  acoustical  wavelength.  The  basic 
dependence  of  the  reflection  coefficient  and  S  upon  <j>  is  sketched  in 
Fig.  Dl.  Large  fluctuations  are  expected.  As  an  alternative  let  us 
consider  the  scattered  and  reflected  sound  as  a  function  of  roughness. 
Figure  D2  illustrates  its  dependence  upon  the  root-mean-square  roughness 
cr  ,  acoustical  wavelength,  and  incident  angle. 


V) 


LARGE  ROUGHNESS 


GRAZING  ANGLE,  ^ 


Fig.  Dl.  Comparison  of  large 
and  small  roughness  as  a 
function  of  <j>  . 


Fig.  Dl.  Comparison  of  large  and  small  roughness  as  a  function  of  <j» 
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Fig.  D2.  For  ffi.  =  1  ,  apparent  reflection  coefficient  ttap  . 

The  conversion  to  apparent  bottom  reflection  is  shown  in  Fig.  D3 
We  assumed  that  the  scattered  coefficient  is  approximately  -20  dB. 


Fig.  D3.  Apparent  bottom  reflection. 
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Often  data  are  expressed  as  bottom  reflection  loss,  in  dB  as  a 
function  of  frequency  and  angle.  Equation  (Bl)  expressed  in  dB  is  as 
follows: 


-10  1og10<R>  =  -[20  1og10  tt]+20[2y2<r2]log10  e 

_2  , 

“  -20  log1Q  HL  +  680  sin  4.  (Dl) 

K 

where  y  =  -£■  sin  <}> 
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